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1. Introduction 

Three-body problem in celestial mechanics is one of pioneer studies on chaotic 
phenomena and, at the same time, one of the oldest and the most famous problems 
in physics. First systematic work has been given by Poincare, for which he won the 
prize established by King Oscar II of Sweden and Norway (Poincare, 1899). Main result 
of his work is that the motion of three particles is very complicated (i.e., what we 
now call "chaos") and then it is generally hard to solve it practically. After Poincare, 
many mathematicians and physicists contributed to this problem (See, for instance, the 
book(Siegel and Moser, 1971) as mathematical literature.). In addition, innovation of 
computers assisted physicists to investigate this problem. However, our understanding 
of three-body problem is far from complete. The investigation is still going on. 

Three-body problem also exists in Coulomb systems. The most popular Coulomb 
three-body system is helium atom which consists of one nucleus and two electrons. In 
early stage of making quantum mechanics, many physicists did not know what to do 
with the spectrum of helium atom, in spite of success of explaining of the spectrum 
of hydrogen atom. For hydrogen atom, stable periodic orbits correspond to eigen 
energies of it. This correspondence is now called the Bohr-Sommerfeld(BS) quantization. 
However, they could not explain the spectrum of helium atom by the BS-quantization 
scheme. Einstein extended the BS-quantization to higher-dimensional classically 
integrable case (now we call it the Einstein-Brioullin-Keller(EBK) quantization.) and 
pointed out indirectly this problem of helium atom (i.e., not integrable case) referring 
to the work by Poincare (Einstein, 1917). About fifty years later, Gutzwiller succeeded 
to derive an approximate quantization condition for hyperbolic cases. His formula 
expresses the density of states in terms of unstable periodic orbits, which is now called 
the Gutzwiller trace formula (Gutzwiller, 1971, 1990). In early 90's, several physicists 
applied the Gutzwiller trace formula to hydrogen negative ion(Gaspard and Rice, 1993) 
and helium atom(Ezra et al., 1991) in collinear configuration, i.e., restricted Coulomb 
three-body problem. Series of these studies showed that these atom and ions in collinear 
eZe configuration are fully chaotic, probably hyperbolic and the Gutzwiller trace formula 
yields nice approximate values for series of eigen energies with angular momentum L = 
numerically(See also the review (Tanner, Richter and Rost, 2000).). Here "e" and "Z" 
represent electron and nucleus, respectively. Thus "eZe" stands for the order of particles 
on a line. Their work has left a question why hydrogen negative ion and helium atom 
in collinear eZe configuration are fully chaotic. We used the term "fully chaotic" in 
the meaning that almost all periodic orbits are hyperbolic. Thus, we used the term 
"fully chaotic" for the system which is almost chaotic, but posseses small tiny island. 
In addition, we shall use the term "hyperbolic" in the meaning that all periodic orbits 
are hyperbolic. 

It was numerically shown that if the mass of nucleus is infinite for eZe collinear 
configuration, the symbolic description of it is complete for the negative energy E < 
and Z > l(Bai, Gu and Yuan, 1998). This result suggested that for the inifinite nucleus 



Classical Coulomb three-body problem 



3 



mass, the dynamics is hyperbolic. 

In this paper, we consider the case that all masses are finite and numerically 
confirm that the dyanmics of two-electron atom and ions in collinear eZe configuration, 
i.e., H , He, Li + , Be 2+ , ... is hyperbolic. Parameters of these systems are the 
mass ratio £ = m n /m e and the charge Z of nucleus, where m n {m e ) are the mass of 
nucleus(electron), respectively. We employ numerical computation for these systems 
with the aid of analytical tool(the triple collision manifold(TCM)) (McGehee, 1974) 
and symbolic dynamics(Tanikawa and Mikkola, 2000) from celestial mechanics. The 
TCM is a manifold which is a set of the initial or final conditions of the triple collision 
orbit, i.e., thus just the triple collision points. The flow on the TCM does not have 
physical reality, because the TCM corresponds just triple collisions. However, the flow 
on the TCM affects the flow near the TCM because of the continuous property of the 
solution for the equations of motion. Thus, we can study the behaviour near the triple 
collision from the flow on the TCM. It is shown that the mass ratio £ plays an important 
role for dynamical behaviour of these systems. In precise, sufficiently large mass ratio 
£ makes the system be fully chaotic. Therefore it is revealed that all of these systems 
are fully chaotic, probably hyperbolic, since the mass ratio is sufficiently large, e.g., 
m p /m e 1840, where m p is the mass of proton. 

We confirm this result as follows: First, we elucidate the ^-dependence of the local 
structure near the triple collision by examining the flow on the TCM. In particular, we 
investigate the stable (or unstable) manifold of two critical points on the TCM. By this 
analysis, we show the change of the triple collision orbit when £ is increased. Let us 
denote the winding number of the stable(or unstable) manifold (on the TCM) of the 
critical point around the body of the TCM by Af. Later we shall define the body of the 
TCM precisely. It is shown that when £ — > , Af becomes infinite and, in the opposite 
limit, when £ is increased, Af is decreased and saturated to certain values. Second, 
we examine the morphology of the Poincare surface of section using the associated 
symbolic dynamics. The symbolic dynamics is naturally introduced by labeling the 
double collision between the particles 1 and 2(2 and 3) by symbol 1(2), respectively. The 
position of the triple collision orbits on the Poincare surface of section is specified. When 
(Z,£) = (1,1), there is a torus on the Poincare surface of section, whose corresponding 
sequence of symbols is . . . 12121212 . . .. The triple collision orbits forms curves on the 
Poincare surface of section. We call these curves the triple collision curves Ctc- in 
this case, Ctc crosses its reversed pair C TC transversely except near the torus. C TC is 
mapped to Ctc by the symmetry operation due to the reverse operation. For helium 
{Z,£) = (l,Am p /m e ), there is no torus. Ctc crosses its pair C TC transversely everywhere. 
From numerical calculation, Ctc{Ctc) ls parallel to the stable manifolds(the unstable 
manifolds) in the Poincare section for the Poincare map, respectively. Therefore, this 
fact is an strong evidence that the stable and unstable manifolds transversely cross each 
other for the case of helium in collinear eZe configuration. In addition, we show the 
relationship between the structure of the stable and unstable manifolds of two critical 
points on the TCM and C T c and C TC on the Poincare section. At this time, it is turned 
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out that the change of the triple collision orbit on the TCM w.r.t £ is related to the 
collapse of the torus whose corresponding symbol sequence is . . . 12121212 . . .. The above 
facts strongly suggest that two-electron atom and ions in collinear eZe configuration is 
fully chaotic, probably hyperbolic. 

The organization of this paper is as follows. In §21 the Hamiltonian of the systems 
considered in this paper is introduced. Due to the attractive interaction between 
neighbouring particles, double collisions surely occur. Therefore, an regularization is 
required before numerical computation. We choose an algorithmic regularization not 
analytical regularization such as Kustaanheimo-Stiefel transformation or Levi-Civita 
transformation. In §H1 the TCM is introduced. The ^-dependence of the flow on the 
TCM is numerically investigated. In particular, the change of the triple collision orbit 
w.r.t. £ is shown. In §H method of symbolic dynamics is applied to the systems for 
the case Z = 1 and £ = 1 and the case of helium. In the former case, there is a torus 
whose symbol sequence is . . . 12121212 . . .. We examine the topological character of this 
system in detail. On the other hand, in the case of helium, there is no torus whose 
corresponding symbol sequence is . . . 12121212 .... In fact, the corresponding orbit for 
the symbol sequence . . . 12121212 ... is an isolated unstable periodic orbit. With the 
result of §3 what this observation means is considered. In £0 the results of this paper 
is summarized. 



2. Hamiltonian and regularization 

We consider three particles 1, 2 and 3 whose masses are mi = m e , m 2 = m n = £m e , 
and m 3 = m e and whose charges are — e, Ze, and — e, respectively. The mass ratio £ is 
given by £ = m n /m e . The Hamiltonian for this system is 

3 v 2 Ze 2 Ze 2 e 2 
H = yJP^_ ^ Ze + ^ (1) 

^2m; |<?i-g 2 | Ifc-ftl Ift-ffsl 
Now we employ the famous scaling for Coulomb systems. 

Qi = ccq'n Pi = Pp'i, E = <yE' (2) 

with 

a = -17773-, P = \E\ 1/2 m e e 2 , 7 = m e e 4 . (3) 

m e e z \E\ 

E is the value of the Hamiltonian H. After this scaling, we get the following 
Hamiltonian. 

_ Pi P'i Pa 1 1 1 

H ~ T + 2£ + Y ~ W^\~ + Z\q[ - q' 3 \ ~ (4) 

For simplicity, we set E' = — 1. Now the particles are arranged in the order q[ < q' 2 < q' 3 . 
If we set (Z, £) = (—1, 1), the system is equivalent to the system which was investigated 
in gravitational three-body problem (Tanikawa and Mikkola, 2000). We set the total 
momentum to be zero and change the variables q'^s to 

Qi = q' 2 - q[, Q2 = <?2> Q3 = <?3 - ?2- ( 5 ) 
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with new momenta 

Pi = -Pi, P2 = p[ + P 2 - P3, A = P' 3 - (6) 

This canonical transformation is generated by the following generating function W. 

W = P 1 (q' 2 - q[) + P 2 q' 2 + P 3 {q' 3 - q' 2 ). (7) 

The final form of the Hamiltonian ri is 

_ Pi Pi P X P 2 11 1 

~ 2^ + 2^ ~ "T" ~Q~i~Q~ 2 + Z{Q 1 + Q 2 ) ~ _1 ' (8) 

with /i = £/(£+ !)■ After the canonical transformation, we replaced (P 3 , Q 3 ) by (P 2 , Q 2 ) 
for convenience. The parameters of this Hamiltonian is the charge Z and the mass ratio 

The potential in Eq.flSJ) has singularities associated to double collisions of 
particles. Thus appropriate regularization is needed. To regularize them, in celestial 
mechanics a transformation such as Kustaanheimo-Stiefel transformation or Levi-Civita 
transformation is usually employed. However, an algorithmic regularization(Mikkola 
and Tanikawa, 1999) is used here. A merit of this choice is that for accurate 
numerical integration of equations of motion, the symplectic integrator method can 
be used(Yoshida, 1990). For usual analytical regularization, the symplectic integrator 
method can not be applied, since there are coupling terms of position and momentum 
in the Hamiltonian, i.e., the Hamiltonian after the transformation mentioned above is 
not a summation of the form ri = T(P) + V(Q). We will use the sixth order symplectic 
integrator method in After regularization, double collision is extended to just elastic 
collision and triple collision is not regularized in general. The interaction between outer 
particles and middle one is attractive, while the interaction between outer particles 
is repulsive. Thus neighbouring particles are always attracted and collide each other. 
Hence the trajectory of three particles is, in general, a sequence of collisions. Among 
collisions, triple collisions exist as very rare events. In fact, the measure of set of 
triple collisions would be zero. However, the triple collisions form the bone structure of 
dynamics of our systems which we will see in and 0] 

Here the symbolic coding of the orbit which will be extensively used in 2] is 
introduced briefly. Let us denote the whole set of orbits in our systems by O. The 
set O is an union of Odc an d Otc-, where One is the set of orbits which only consist 
of double collisions and Otc is the set of orbits which include triple collisions. 

O = O dc U O tc . (9) 

For an orbit of Odc, it is natural that the double collision with the particle 1 and 2(2 
and 3) is labeled by symbol 1(2), respectively. Therefore, the symbol set for Odc is 
A = {1, 2}. Let us express an given orbit of Odc as a sequence of symbols as follows: 

n = . . .n_ 2 n_i • n nin 2 . . . , (10) 

where n, G A ,i e Z. In order to describe an orbit of Otc as a symbol sequence, we 
need another symbol, i.e., 0. We set A = {0, 1,2}. The triple collision orbit is started 



Classical Coulomb three-body problem 



6 



and/or ended by the triple collision. Let us label the triple collision by symbol 0. The 
orbit which experience the triple collision can not be continued. Thus we regard as 
continuing endlessly after the first in the future and before the last in the past, if 
appears in the sequence. For example, a triple collision orbit is represented as 

n = . . . 000n_ / n_ (/ _ 1 ) . . . n_i ■ n nin 2 (11) 

or 

n = . . . n_^n_2 n -i • n o n i ■ ■ ■ ^mOOO . . . , (12) 

or 

n = . . . 000n_ ; n_(/_!) . . . n_i ■ n ni . . . n m 000 . . . , (13) 
where G A . The shift operator a on A z is defined by 

<x(. . . n_ 2 n-i ■ n nin 2 ...) = ... n_ 2 n_ l n Q ■ n x n 2 (14) 

We sometimes call finite symbol sequence the word. For instance, 12 and 1211 are words. 



3. Triple collision manifold 

The triple collisions are, in general, essential singularities and thus are not regularized. 
They are rare events. Therefore, it is hard to visualize the triple collisions. In order 
to investigate the structure near the triple collision, technical method is required. For 
celestial problem, McGehee has developed such method(McGehee, 1974). He has derived 
the equations of motion for the flow just on the triple collision. Its derivation is successive 
application of tricky transformations to the equations of motion Eq. (jlj) and the energy 
conservation relation Ti = T(p) — U (q) = E. A manifold on which the orbit experiences 
just the triple collision is obtained by setting the moment of inertia to be zero (i.e., 
just triple collision) in the final energy conservation relation. This substitution is 
meaningful owing to tricky transformations (i.e., scalings and time-transformations). 
This manifold is called the triple collision manifold(TCM) . Thanks to similarity between 
celestial problem and Coulomb problem, this method can be also applied to our Coulomb 
problem. Since the transformations are tricky and complicated, in order to be self- 
contained, we show the derivation of it rapidly in the Appendix. For the detail of the 
derivation, the readers are recommended to consult with the article (McGehee, 1974). 
Starting from the energy conservation Eq.(J3J and the equations of motion for Eq. (|4*j) . 
after lengthy calculation( i.e., six times of changes of variables ), we obtain the energy 
conservation relation with r = 0, i.e., the moment of inertia r 2 = q\ + £g| + q\ is zero 
(just the triple collision) or with alternately the total energy E = f : 

w 2 + s 2 -l + (l-s 2 ) 2 VV(s)-V = 0. (15) 

Equation (JT5|l defines a surface in (s, v , ^-coordinates, i.e., the TCM. The TCM is 
topologically equivalent to a sphere with four holes. Schematic picture for the TCM 

% r and E appears as a term rE in the energy relation and equations of motion(See the Appendix.). 
Therefore, the TCM and the flow on the TCM are the same for either r = and E = 0. 



Classical Coulomb three-body problem 



7 



is depicted in Fig. [TJ For the definition of s,v,w,W(s), see the Appendix. We only 
explain what s represents, s represents the configuration of three particles, s is valued 
in [—1, 1]. s = — l(s = 1) corresponds to the double collision between the particles 1 and 
2(2 and 3), respectively. The flow on the TCM is determined by the following equations 
of motion 



ds 
dr 



1-s 2 
W(s) 



dw 2s(l - s 2 ) „ 1 W'(s) , , „. A 1 - s 2 

For the definitions of X,s,v,w, and r, see the Appendix. There are two critical points 
c = (0, — v c , 0) and d = (0, v c , 0), two infinite arms(a and b), and two infinite legs(e and 
/), where v c = W(0)^. Most of orbits on the TCM comes from infinity of one of two 
legs winding around it and goes out to infinity of one of two arms winding around it. 

There is an important property for the orbits with r ^ (i.e., for total flow). 

Property 3.1: For TC = E < 0, the orbits runs inside of the TCM, i.e., 
w 2 + s 2 + (1 - s 2 ) 2 W(s)~ 1 v 2 < 1. On the other hand, for H = E > 0, the orbits 
runs outside of the TCM, i.e., w 2 + s 2 + (1 - s 2 ) 2 W(s)~V > 1. 

Here after we consider the case of 7i = E < 0. 

If the orbit crosses s = — l(s = 1), the corresponding trajectory in the configuration 
space experiences the double collision between the particles 1 and 2(2 and 3) with the 
symbol 1(2), respectively. After some collisions they are going out to infinity by changing 
its binary or not. There also exist the triple collision orbits. By definition, the triple 
collision orbits are the orbits which start/end at the TCM(i.e., r = 0). This means that 
the triple collision orbit can not be regularized, that is, the triple collision orbits can 
not continue after/before the triple collisions. As shown by McGehee(McGehee, 1974), 
the triple collision orbits form one parameter family. This fact is understood from the 
stability analysis of the critical point c and d. The critical points c and d are the fixed 
points of the flow Eq. pfij) . At the same time, they are the fixed points of the total flow 
Eq.([46|) which is not restricted to r = 0. The stability analysis of the fixed points c and 
d shows that dim(W s (c)) = 2,dim(W"(c)) = 1 and dim(W s (rf)) = 1, dim(W M (rf)) = 2, 
where W s (x) and W u (x) are the stable and unstable manifolds of x, respectively. In 
FigEl the schematic picture of W s (c), W"(c), W s (d), and W u (d) is depicted. One branch 
of W s (c)(we call it >V s '*(c)) comes into c on the TCM from the outside of the TCM(i.e., 
r ^ 0) along the f-axis. Similarly, one branch of W u (<i)(we call it W u '*(d)) goes out 
d on the TCM to the outside of the TCM(i.e., r ^ 0) along f-axis. The other branch 
of W s (c) and W(c) (we call them Wtcm{c)) runs on the TCM and winds around 
the TCM. Similarly, the other branch of W a (d) and W s (d) (we call them W T CM(d)) 
runs on the TCM and winds around the TCM. In the outside of the TCM(i.e., r ^ 0) 
near the critical point c and d, there exists one parameter family of the orbits which 
approach to Wtcm(c) as r — ► oo. Similarly, there exists one parameter family of the 
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orbits which escape from WxcM^d) as r — > oo. Therefore, Wtcm{.c) and Wtcm(^) 
determine the behaviour of the triple collision orbits. Unfortunately, this discussion 
here is limited to the neighbourhood of the critical points c and cf(i.e., local property). 
In the next section, we investigate Wtcm(c) and Wtcm(^) numerically to show the 
global topological property of the triple collision orbits. For later use, we call the 
part of the TCM between two critical points c and d the body of the TCM, in precise, 
{(s, v, to); — v c < v < v c , (s, v, w) on the TCM}. 

Our interest in this section is focused on the ^-dependence of the flow on the TCM. 
In the celestial problem in collinear configuration considered in the section 10 of the 
article(McGehee, 1974), it is shown that for the case of m± = m^ = m and m 2 = em, 
when e — > 0, the orbits on the TCM wind around the TCM infinitely. Thanks to the 
similarity between celestial problem and Coulomb problem, for our Coulomb systems, 
the same argument is easily shown following the discussion of the article(McGehee, 
1974). 

Property 3.2:When £ — > 0, the orbits on the TCM wind around the body of the 
TCM infinitely often. 

This is due to the fact that 4^ — > as £ — > 0. The proof is the same as the 
Proposition 10.1 of the article(McGehee, 1974). 

In order to show the ^-dependence of the flow on the TCM for large £, we 
numerically calculate Wtcm(c) and Wtcm(^)- For this numerical integration of Eq. ([TKjl . 
we used the fourth order Runge-Kutta method. Figure 01 depicts Wtcm{.c) and Wtcm(^) 
for (Z,£) = (1,1). Figures Ufa) and (b) depict one branch(the unstable manifold) of 
Wtcm(c) for the cases (Z,£) = (1, 0.1) and (1, 6), respectively. The change of Wtcm{c) 
w.r.t. £ is clearly seen in FigEJ When £ is increased, the winding number M of 
Wtcm( c ) or y^TCM^d) around the body of the TCM is decreased. Although further 
£ is increased from £ = 6, further change does not occur. Thus, when £ is increasing, 
M is monotonically decreasing to a certain value. This property is also observed for 
different values of Z. 

Remark: For our Coulomb system, the potential part of the Hamiltonian does 
not include the masses. Thus, v c is independent of £. For the gravitational case, the 
TCM has the same topology compared with our Coulomb system. But the potential 
part of the Hamiltonian includes the masses. This makes difference in the behaviour of 
the flow on the TCM. If m\ = l,m 2 = £, and m 3 = 1, v c ~ as £ — > oo. Therefore, 
when £ — > oo, the winding number of Wtcm(c) or Wtcm(gO may not decrease as well 
as in the case of our Coulomb system. Thus, it is supposed that in the gravitational 
three-body problem in collinear configuration, the system hardly becomes hyperbolic, 
namely it always has stable orbits(i.e., tori). 

To summarize, when £ — > 0, the orbit on the TCM wind around the body of the 
TCM infinitely often. When £ is increased, the number of times that the orbit winds 
around the body of TCM is monotonically decreasing to a certain value. The latter fact 
is related to what we examine in the next section. 
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4. Symbolic dynamics 

In this section, we examine the global structure of the Poincare surface of section by 
using the associate symbolic dynamics. The Poincare surface of section is defined as 
follows: We denote the position on the line Q\ = Q 2 in the (Qi, Q 2 )-pl&ne by R, i.e., 
R = Qi = Q 2 . On this line, the momenta P\ and P2 are specified by new variables 9 
and R. 



Pi = Vfcos9 + J J— -sinO, 

n + 2 

P 2 = VTcos6 - ] J^L S m6, (17) 

where T is the kinetic part of the total energy. Our surface of section is the (9, i?)-plane 
where < R < R ma x, < 9 < 2n and R ma x = 2 — j^. We denote this plane by T>. We 
define the map x from the point z G T> to the bi-infinite symbol sequence n which is an 
itinerary of the orbit started from z at the time zero. 

X (18) 
z 1— > n = . . . n_ 2 n_i ■ n nin 2 . . . 

We define the map y/ + ) from the point z to the semi-infinite symbol sequence n^. 

y(+) ■ J) -> A N 

' <+) 19 

z 1— > n w = n_i • n nin 2 . . . 

We also consider the following map from the semi-infinite symbol sequence to the 
finite symbol sequence with length / + 2: 

S;(n_i • n ni . . .) = n_i • n ni . . . n/. (20) 

Let us introduce 

V 1 = {(9, R) : < 6 < 71, < R < R max }, 

V 2 = {(9,R):n<9<27t,0<R< R max }. (21) 

There are two symmetry operations on T>. (1) if the orbit starting z = (9,R) G T>\ 
has the symbol sequence . . . m^ 2 m^\ ■ n§n\ . . ., then . . . iiiUq ■ m_im_2 ... is a symbol 
sequence for the point (2ir — 9,R) e T> 2 . (2) if -nin 2 n^ . . . is the future symbol 
sequence corresponding to the orbit starting at {9, R) G T>\ and if •mim 2 m^ ... is 
the future symbol sequence corresponding to the orbit starting at (ir — 9,R) G T>\, 
then . . . m'^m^mi ■ nin 2 n 3 ... is the bi-infinite sequence corresponding to the full orbit 
starting at the point (9, R) G T> x , where m! i = \\im i = 2 and — 2 if rrii — 1. Thanks 
to these two symmetry operations, in order to investigate the global structure of V, it 
is sufficient to study only the future orbits for the points in T>\. 

We numerically construct the map x^ i n the following way: we consider the 
rectangle lattice whose lattice size is 1800 x 1000 for T>\. For each lattice point 
z nm, 1 < n < 1800, 1 < m < 1000, we numerically obtained the truncated symbol 
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sequence S 60 (x( + )(z nm )) by integrating the equations of motion for the Hamiltonian 
Eq.fjHJ) using the sixth order symplectic integrator method as mentioned in ^21 

As established for celestial problem(Tanikawa and Mikkola, 2000), the following 
properties are also true for our Coulomb system: 

Property 4.1: A trajectory in the (Qi, Q2)-plane transversely crosses the line 
Qi — Q2 except at (9, R) = (0, 0), if it does at all. 

Property 4.2: If a trajectory crosses the line Qi = Q2 on the (Qi, Q2)-plane, a 
double collision occurs before the trajectory again crosses it. 

The proof for the above Properties is the same as in the article(Tanikawa and 
Mikkola, 2000). 

From the latter Property 4.2, if we put the initial condition z in T>i whose trajectory 
crosses the line Q\ = Q2, then the corresponding orbit has the symbol sequence 
. . . ?2_3n_2l • 2niU2 . . ., where n« G A. 

The following two things are numerically checked for our Coulomb systems (i.e., 
the case of (Z,£) = (1,1) and the case of the helium) as observed for celestial 
problem(Tanikawa and Mikkola, 2000): The plane T>i is divided into two regions of 
the points z having Si(%( + )(z)) = 1-22 and 1 • 21. Furthermore, it divided into four 
regions of the points z having S 2 (x (+) ( z )) = 1 • 222 > 1 • 221 > 1 • 212 and 1 ' 211. This 
procedure is repeatedly applied. Then T>\ is divided into smaller regions. The second 
thing is that the boundaries of regions of different symbol sequences form curves in T>\. 
It is turned out that these curves are initial conditions of orbits which end in triple 
collision. So we call these curves the triple collision curves. We denote them by Ctc- 

Figures El and IH1 depict the triple collision curves Ctc f° r the case of (Z,£) = (1, 1) 
and the case of the helium, respectively. For the case of (Z, £) = (1, l)(FigJH|), there 
is a torus in T>\. The simple stable orbits(the torus) may correspond to the Schubart 
orbits in celestial problem(Schubart, 1956). Therefore, we call these stable orbits the 
S-orbits. In the context of atomic physics, it is usually called the asymmetric stretch 
orbit which was found by Simonovic and Rost(Simonovic and Rost, 2001). In Fig|3 
the region of the torus(S-orbits) is shown as a triangle area located at 9 — f . For 
the point z in this triangle area, it is that x^K 2 ) — 1 ■ (21)°°. In Fig|Sl we show the 
regions of the points z whose symbol sequence is H 6 (%( + )(z)) = 1 ■ 2nin 2 . . . n 6 , where 
Hi G Aq(z = 1, 2, . . . , 6). There are the missing regions of the points z whose symbol 
sequence S;(x^ + ^(z)) includes the words 1122 and 2211. Our numerical calculation up 
to the word length 15 showed that the inadmissible words are only the words including 
the words 1122 and 2211. Note that 1122 and 2211 are also inadmissible words for 
the celestial problem(the case of (Z,£) = (—1, 1)) and other inadmissible words exist 
(Tanikawa and Mikkola, 2000). On the other hand, for the case of the helium(FigE|), 
i.e., large £, there is no torus. Figure |H1 shows the regions of the points z whose symbol 
sequence is H 6 (x^ + ^(z)) = 1 • 2n\ri2 ■ ■ ■ n e , where rij G A (i = 1, 2, . . . , 6). In this case, 
there is no missing regions(i.e., no inadmissible words). When the symbol sequences 
1.2ni . . . rii with length I + 2 are considered, T>\ is divided into 2 l partitions. 
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In order to examine the hyperbolicity of the system, we investigate the foliated 
structure of V>\. Using the second symmetry of T>i, we construct the triple collision 
curves whose orbit is started at triple collision in past. We denote the triple collision 
curves obtained in this way by C TC . In FigsJ7|and|Bl we depict the triple collision curves 
whose orbit is started and/or ended at triple collision for the case of (Z,£) = (1, 1) 
and the case of the helium, respectively. Figsd and |H1 are constructed from FigsEl and 
El It is clearly seen that for Fig|H] Ctc an d C TC transversely cross each other, while 
for Fig0 Ctc and C TC transversely cross each other except near the torus. Since the 
dynamics of our system is continuous, we expect that Ctc an d C TC do not cross the 
stable and unstable manifolds. For discontinuous system, such as billiard system, this 
is not the case. In fact, it is numerically confirmed that when £ is sufficiently large, 
the triple collision curves Ctc{Ctc) * s parallel to the stable(unstable) manifolds in the 
Poincare section for the Poincare map, respectively. Therefore, FiglH^sufficiently large 
£) manifests that the dynamics of the helium in collinear eZe configuration is hyperbolic. 
With some parameter values when the torus exists(for small value of £), Ctc and C TC 
do not foliate. In this case, it is observed that the tangency of the Ctc an d C TC . This 
may manifest the tangency of the stable and unstable manifolds. 

The relationship between the observation in the previous section and the 
observation in Figures [7| and |H] is unclear at present. Next we clarify this relation by 
transforming the Poincare plot in (8, i?)-coordinates into that in (s, v, u>)-coordinates. 
The Poincare surface section in (s, v , w )-coordinates corresponding to that in (8, R)- 
coordinates is just the plane s = 0. The transformation from (8, R) to (v, w) is as 
follows: 

1 



w 



*V2 
l V2 



2V2 ( 2 




(22) 



where 



P[ = cos 8 



e + 2 



sin 8, 



P' 



and 



cos 8 — 

( -Pi 
Pi 



\ 



■ sin 6, 



P 1 



\ 



Pi 



J 



1 

72 





V 1 J 



(23) 



(24) 



Figure El shows the Poincare plot of the triple collision orbits in (v, u>)-coordinates, 
which corresponds to Figure [3 in (8, i?)-coordinates for (Z, £) = (1,1). w > 0(w < 0) 
corresponds to U 1 (V 2 ), respectively. First, we note that there is a special solution along 
the line s = w = 0, which starts from the critical point d and ends in c. It is given by 



f r tanh 



TO 
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r = 




(25) 



where r is determined by the initial condition. This solution shows that W s '*(c) and 
W u, *(d) are degenerated. This type of solution also appears in symmetric gravitational 
four-body problem(Sekiguchi and Tanikawa, 2002). Second remarkable point is that 
the triple collision curves Ctc and its reversed pair C TC accumulate at ten points on 
the TCM(i.e., w 2 + W{fi)~ x v % = 1). Figure EH for = (1,1) shows one branch 

of Wtcm(c). As shown in FigUHl these points are the points at which Wtcm(c) and 
y^TChi{.d) cross the plane s — 0. We denote these points by Vtcm,t=q- It is easily 
understood that the number of points of Vtcm,t=o is related to the existence of the tori in 
the Poincare surface of section s = 0. If the tori exists, its outer most torus has periodic 
points. These periodic points have the stable and unstable manifolds. Some branches 
of these stable and unstable manifolds run toward Wtcm(c) and Wtcm(^)- Therefore, 
the number of the points of VrcM,r=o is directly related to the period of the periodic 
points associated to the tori. At the same time, the number of the points of P TCM ,r=o 
just corresponds to the winding number Af of Wtcm( c ) or VVtcm(^) around the body 
of the TCM as observed in the previous section. In FigHH the case of (Z,£) = (1,7) 
is shown. As the result of the previous section, when £ is large enough, the winding 
number of Wtcm(c) and WrcM^d) saturates to certain value. In other words, when £ is 
large enough, the number of the points of Vtcm,t=o also saturates to certain value which 
is, in fact, six(i.e., no torus case). Thus, the existence of torus in the Poincare section 
s = is monitored by the number of points of Vtcm,t=q- When £ is large enough, the 
triple collision curves Ctc transversely cross C TC . This would be a strong evidence of 
hyperbolicity of the system with large £, since it is numerically confirmed that Ctc(Ctc) 
is parallel to the stable(unstable) manifold for the Poincare map, respectively. 

The critical value of £ C {Z) at which the winding number M is minimized, is 
calculateted. In Tabled we summarize the result. For Z = 1,2,3,4,5, the critical 
value £ C (Z) is order of O(10). These critical value is numerically obtained. Thus, 
these are not true critical values. However, we expect that at these true critical value 
the unstable manifold of c and the stable manifold of d degenerate. At these critical 
value £ C {Z), the tori disappear. At present, we do not know precise mechanism of this 
disapearance of the tori, e.g., even whether KAM-scenario is applied or not. 

From the above numerical observation, we can state a conjecture. For the system 
Eq-dU) with (Z,£), there exists the critical value £ C (Z) such that for £ > £ C (Z) 
the system is hyperbolic. 

Since for actual two-electron atom or ions the mass ratio is large, i.e., m p /m e w 
1840, this means that the classical dynamics of H~,He,Li + ,Be 2+ ,. . . in collinear eZe 
configuration is hyperbolic. For inifinite mass ratio, it was numerically shown that the 
symbolic description is complete and the dynamics probably hyperbolic(Bai, Gu and 
Yuan, 1998). 
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4.1. The case of{Z,£) = (1,1) 

In this subsection, we investigate the detailed structure of the Poincare surface of 
section for the case of (Z,£) = (1,1). In this case, there is a torus whose symbol 
sequence is 1 • 21212 ... as shown in FigsJHl and [7| This torus has the periodic points 
with period 6 as the outermost part. Note that we include the periodic points in T>2 
and count the period. These periodic points a, (3, 7 has stable and unstable manifolds. 
a : (6, R) = (0.5?r, 1.385), (5 : (0, R) = (0.5292tt, 1.275),and 7 : (9, R) = (0.4708vr, 1.275). 
From the numerical calculation, one branch of W s (a) is equal to one branch of W n (/3). 
In the same way, one branch of W s (j3) is equal to one branch of W n (7) and one branch 
of W s (7) is equal to one branch of W u (a). We can construct the stable manifolds of 
a, (3, and 7 outside the S-orbits by examining the long orbits with the symbol sequence 
1-21212 • • •. The result is depicted in Fig lT2T a). In Fig IT2T b). we also depict the unstable 
manifolds by using the second symmetry of T>\. The stable manifolds of a, (3, and 7 
outside S-orbits are basic boundaries of the partitions of T>\. 

We further examined the symbol sequence S 60 (x ( ' + ' ) ( z mn.)) for each z mn (0 < 
m < 1800,0 < n < 1000). As a result, V 1 is divided into, at least, ten partitions: 
SJiii = 1,2, ...,9). S stands for the S-orbits. Other partitions = 1,2, ...,9) 
are defined as follows. We examined the symbol sequences S 60 (x^ + )(z)) along the line 
9 = 0.27T, 0.57T, 0.557T and 0.77T. For the line 9 = 0.2ir, we find that the symbol sequences 
are distributed as it decreases from 1 • (2)°° to 1 ■ (21)°° with increasing R. We divide 
symbol sequences into three groups: 1 • (2)°°, 1 • (2) n 12 . . ., and 1 • (21) n . . . (n > 2). We 
call the regions with these symbol sequences J7, and I\. 

For the line 9 = 0.5tc, the symbol sequences are distributed as it increases from 
1 ■ 2(1)°° to 1 ■ (21)°°. We divide symbol sequences into three groups: 1 ■ 2(1)°°, 
1 ■ 2(l) n 21 . . ., and 1 • (21)" . . .(n > 2). We call the regions with these symbols sequences 
I 8 , h and I 2 - 

By examining the symbol sequences along the line 9 = 0.557T and 0.77T, we find that 
along the line R = 1.3 from 9 = tt to 0.57T, the symbol sequences are distributed as it 
decreases from 1 ■ 21(2)°° to 1 • (21)°° with decreasing 9. We divide symbol sequences 
into three groups: 1 ■ 21(2)°°, 1 ■ 21(2) n 12 ... and 1 • (21) n . . .(n > 2). We call the regions 
with these symbols sequences Ig, Iq and I3. FiglT37a) shows these partitions of T>\. In 
Fig lTBT b). the time- reversed partitions by the second symmetry of T>\ are shown with 
them of Fig ll3r a). In Table the characteristic feature of 5* and liii = 1, 2, . . . , 9) is 
summarized. 

Next we prepared many points in each partition and checked which partitions they 
are mapped to. In Table El we summarized the transition among the partitions. From 
Table El we drew the diagram of the transitions among partitions in FiglTH From this 
diagram, we know that typical orbits travel around the regions near the torus(Ji,/2, and 
^3) and/or around the regions in which the orbits feel large instability (145/5 an d ^6)) 
and, that in most cases, they escape to the regions(/ 7 , 1% and J 9 ) except the S-orbits 
and non- wandering orbits(probably repellor). 
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4-2. The case of helium 

For the case of helium, Figs. El and |H1 tell all features of the partitions. The partitions 
formed by Ctc is, at first sight, expected to be the Markov partitions. If we obtain 
Ctc using the symbol sequences 1.2n\ . . .ni, then D\ is divided into 2 l partitions. From 
numerical observation, there is no inadmissible words. However, it is clear that 1 and 
2 are removed, since these orbits correspond to fixed points at inifinity. The Markov 
partition should be shaped as parallelepiped. Thus probably, in mathematical rigorous 
sense, the partitions in our Poincare section are not the Markov partitions. If we 
construct the Markov partitions in our Poincare section, they would be infinite. This 
reflects the fact that the dynamics of helium atom in collinear eZe configuration exhibits 
intermittency (Richter, Tanner and Wintgen, 1993, Tanner and Wintgen, 1995, Tanner, 
Richter and Rost, 2000). The intermittency observed in helium atom in collinear eZe 
configuration is the behaviour that the electrons are going back and forth between long 
flight (i.e., almost ionization) and short flight (i.e., successive collisions with the nucleus). 
This intermittency is due to the existence of the fixed points at infinity. Its behaviour 
is very similar to hyperbola billiard(Sieber and Steiner, 1990). The transverse crossing 
of Ctc and C TC strongly suggests that the dynamics is hyperbolic. 

4.3. The Poincare plots for other two- electron ions 

For the case of the inifinite nucleus mass, it was investigated by Bai et al(Bai et al, 
1998). It was numerically shown that for Z > 1, the symbolic description is complete 
and it suggests that the system is hyperbolic. We now consider the finite mass case. The 
investigation above does not consider the Z-dependence of the dynamics of two-electron 
atom and ions in collinear eZe configuration. Now in order to strengthen our claim that 
the dynamics of them is fully chaotic, probably hyperbolic, we calculated the Poincare 
plots for thirty atom and ions among the systems H~, He, Li + , ... , Fm 98+ in collinear 
eZe configuration. We do not plot them since there is no space. This calculation showed 
that even though Z becomes large at the order 0(100) and then the interaction terms 
becomes small in Eq.fjHJ), their large mass ratio £ overcomes the Z-dependence to result 
in the fact that the Poincare plot is filled by chaotic sea except the escape region and 
there is no visible torus. These observations strengthen our claim. 

5. Experimental aspect 

From our observations, we know that the system with small mass ratio qualitatively 
differs from the system with large mass ratio, namely the existence of stable orbits. 
This difference would be experimentally observed. A possible candidates are the 
antiproton-proton-antiproton (p-p-p) system, the positronium negative ion(Pr-(e-e-e)). 
which corresponds to the case of Z = 1,£ = 1. We neglect relativistic effects, 
bremsstrahlung and hyperfine interaction. If this system has the bound states, the 
eigenenergies possesses the effects of the torus. Most convenient analysis is the Fourier 
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transform of the density of states for the spectrum with the angular momentum L = 0, 
which gives the information on the length and the stability of periodic orbits. For 
the positronium negative ion, the EBK quantization was done(Simonovic and Rost, 
2001). Stable antisymmetric orbit was obtained and was quantized to explain the energy 
spectrum. From our observation, as a ZeZ configuration, the positive hydrogen molcule 
H2 have stable periodic orbits. 

In addition, although large Z two-electron ions are experimentally unrealistic, 
small Z two-electron ions serve us the data for the manifestation of our finding (i.e., 
hyperbolicity (strong chaotic property)). and He have been already analyzed by 
the semiclassical method (Gaspard and Rice, 1993, Ezra et al., 1991). Therefore, Li + is 
another candidate. If we change Z, the degree of the intermitency of the dynamics would 
change(Tanner and Wintgen, 1995). This change may be reflected in the behaviour of 
the quantum defect of H~, He, Li + , Be 2+ , .... 

6. Summary 

In this paper, we have investigated the bifurcation in classical Coulomb three-body 
problem in collinear eZe configuration with finite mass of all three particles. In 
particular, the main result is that when the mass ratio £ is changed, the change of 
the flow on the TCM is directly related to the existence (or collapse) of the tori. This 
result suggests that if the mass ratio £ is sufficiently large, the dynamics of these system 
is hyperbolic. This result is consistent with the result (i.e., infinite nucleus mass) by Bai 
et al(Bai, Gu and Yuan, 1998). But by our analysis, it was shown that there surely 
exists a threshold value of £ for hyperbolicity. These threshold value is order of O(10). 
Therefore, H~, He, Li + , Be 2+ , ... in the collinear eZe configuration is hyperbolic. Since 
our analysis is based on numerical one, we do not have the proof of this fact. We call 
for rigorous proof of this result. 
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Appendix 

In this Appendix, we derive the TCM for our Coulomb systems. Its derivation is 
completely parallel to the case of celestial problem. The readers are recommended 
to consult with the article for detail (McGehee, 1974). We start from the Hamiltonian 
Eq.®. 

H = T(p) - U(q) = E. (26) 
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where we removed the dash in Eq.(J3J). T(p) is the kinetic part of 7i and — Z7(q) is the 
potential part of 7i. The equations of motion in the Cartesian coordinates are given by 



dp 

~dt 



W(q), 



with the mass matrix 

" 1 
M= £ 
1 

We set the center of mass to be zero. 



(27) 



(28) 



(29) 



Ql + &2 + <?3 = 0. 

We set the total momentum to be also zero. 

Pi + P2 + Ps = 0. (30) 

We consider the transformation from the Cartesian coordinates (qi, q 2 , qziPiiP^iPz) with 
Eqs. ()29|) and (J3UJ) to the McGehee's variables (r,s,v,w). First, we set 

r = (q T Mq)^ 2 , 



s = r q. 



(31) 



r 2 is the moment of inertia, s represents the configuration of three particles. The triple 
collision corresponds to r = 0. Now we parameterize the variable s by single variable s. 
To do so, we set 



where 



a = (ai, a 2 , a 3 ), a\ = a 2 < a 3 
b = (61,62,63), 61 < 6 2 = 63 

a T Ma = b T Mb = 1. 



(32) 



(33) 



Since we fix the center of mass, namely a± + £0,2 + 03 = and 61+^62 + 63 = 0, then we 
have 



1 



1 



1+e 



(1 + 0(2 + 



1 + 0(2 + y/(l + 0(2 + 0, 
1 1 



(34) 

;i + 0(2 + ' v / ( 1 + 0(2 + ' v / ( 1 + 0(2 + 0: 

It is shown that the variable s is parametrized by the variable s G [—1, 1] as follows (See 
the article for detail(McGehee, 1974).): 

s = (sin(2A))" 1 [sin(A(l - s)a + sin(A(l + s))b] , (35) 

with 



cos(2A) = a T Mb 



1 



(36) 
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Here A is valued as < A < |. We denote the map from s G [— 1, 1] to the configuration 
s by S: 

S: [-1,1] - {s:s T Afs = l, Sl + ^ 2 + S3 = 0}, 
s i— > s. 

Note that s = — l(s = 1) corresponds to the collision of the particles 1 and 2(2 and 3), 
respectively. U(s) becomes 

V(s) = U(S(s)) 

r i i 

~ Sm(2 } [(62-6i)sin(A(l + s)) + (a 3 -a 2 )sin(A(l-s)) 

1 

~Z{(b 2 - h) sin(A(l + s)) + (a 3 - a 2 ) sin(A(l - s))} 
Here we set two matrices Ai and A 2 



(38) 





1 


1 


1 " 







1 


-1 






1 


1 


1 


, A 2 = 


-1 





1 


(39) 




1 


1 


1 




1 


-1 








and define the matrix A 

1 



A 



2 + ^ 



A X M + 



M^Ao. 



(40) 



Now the variable s is defined as the inverse of Eq. (J35|) . 

^^(s). (41) 
We further set 

I T 

v = r 2 p s, 

w = r l(l - s 2 )W(s)-h T A T p, (42) 

where 

W(s) = 2(1 - s 2 )V(s). (43) 

We also employ time-transformation two times. 

dt = r§dt' and dt' = A(l - s 2 )W(s)-^dT. (44) 

In short, the change of variables Eqs.(jHH),(jJTJ),(|I2J), and (|4^jl give the wanted relations, 
i.e., the energy conservation and the equations of motion. The energy conservation 
becomes 



1 - 



2w 2 



2(1 -s 2 



1 - s 2 W(s) 
The equations of motion become 
dr A(l - s 2 ^ 



{v 2 - 2rE) - 1. 



(45) 



dr ~~ W(s) 1 / 2 
dr 2 K ' 



H/(s) 



(w 2 - ArE) 
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ds 

dr 
dw 

dr 



2s(l-s 2 ) 
S+ W(s) 

Ml"* 2 ) 



(v 2 - ArE) + 



1 W'{s) 

2 W(s) 



(l-s 2 -w 2 ) 




(46) 



Finally, we set r = 0. Thanks to the above transformations, this substitution is 
meaningful, since the singularities from the double collisions have been removed. From 
Eqs. fll3|) and (JlT)|) . we obtain Eqs.(JTSJ) and (|TB|). We can also set E = 0. Then the 
same energy relation and the equations of motion are obatained, since r and E appear 
as a term rE. The dynamis of scattering flow for E — was investigated by Bai et 
al(Bai,Gu,Yuan, 1998) for the case of the inifmite nucleus mass. 
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Figure 1. The triple collision manifold: There are two critical points c and d. 



Figure 2. The schematic picture of the stable and unstable manifolds of the critical 
points c and d: W(c) and W(d) are depicted. dim(W s (c)) = 2, dim(W u (c)) = 1, 
dim(W s (d)) = 1, dim(W u (d)) = 2. 



Figure 3. The stable and unstable manifolds of the critical points c and d on the 
TCM for (Z,£) = (1, 1): Wtcm(c) and WrcM{d)- Two circles indicate the positions 
of two critical points c and d. 

(a) (b) 



Figure 4. The unstable manifold of the critical point c on the TCM: (a)(Z, £) = 
(1,0.1), (b)(Z,0 = (1,6). 



Figure 5. The triple collision curves Ctc f° r the case of (Z, £) = (1, 1): We also plot 
the torus region (triangle area). The curves are obtained by the symbol sequences with 
length 7 in future. 



Figure 6. The triple collision curves Ctc f° r the case of helium: The curves are 
obtained by the symbol sequences with length 7 in future. 



Figure 7. The triple collision curves Ctc and C TC for the case of (Z, £) = (1, 1): We 
also plot the torus region(triangle area). These curves are obtained from the data of 
Fig.5. 



Figure 8. The triple collision curves Ctc an d C TC for the case of the helium: These 
curves are obtained from the data of Fig. 6. 



Figure 9. The Poincare plot of the triple collision orbits Ctc an d C TC in (v,w)- 
coordinates for (Z, £) = (1, 1). It corresponds to Fig. 7. 



Figure 10. The Poincare plot of the triple collision orbits and one branch of 
the unstable manifold of the critical point c for (Z,£) = (1,1). Three crossing 
points(including the critical point c) where the unstable manifold crosses the Poincare 
section s = are indicated by squares. 
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Figure 11. The Poincare plot of the triple collision orbits Ctc and C^c f° r 
(Z, £) = (1,7) in (v, w)-coordinates. 



(a) (b) 



Figure 12. The stable and unstable manifolds of a, /3, and 7 for the case of (Z,£) = 
(1,1): (a) The stable manifolds of a, (3, and 7, (b) The unstable manifolds (dotted 
lines) of a, (3 and 7 are also added to (a). 



(a) 
(b) 



Figure 13. The partitions in D\ for the case of (Z, £) = (1, 1): (a) The partitions in 
T>i, (b) The time-reversed partitions are also added to (a). 



Figure 14. The transition diagram for the case of (Z, £) = (1, 1) among the partitions 
derived from Table 2. 
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1 


6.4 


2 


15.0 


3 


23.5 


4 


32.0 


5 


40.4 



Table 1. The critical value £ C (Z). 



Partition 


Character 


Symbol sequence 


S 


S-orbits(torus) 


1 • (21)°° 


h 




1- (21) n ...(n> 2) 


h 




1- (21) n ...(n> 2) 


h 




1- (21) n ...(n> 2) 


h 




1 • (2) n 12...(n > 2) 


h 




1 • 2(l) n 21 . . . (n > 2) 


h 




1 • 21(2) n 12. . . (n > 2) 


h 


Escape 


1 • (2)°° 


h 


Escape 


1 • 2(1)°° 


h 


It is mapped to h- 


1 • 21(2)°° 



Table 2. The characteristic feature of each partition for {Z,£) = (1,1). 



Transitions 
S -> S 

h ~» -^2, hi h 

h ~» -^3) -^6; -^9 

/ 3 ^h 

h —> hi h, hi hi hi h 
h ~> hi hi hi hi hi h 

h^h 

h — > Escape 
Is — ► Escape 
/ 9 ^/r 



Table 3. The transitions among the partitions for {Z,£) = (1,1). 



